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Abstract. Recovering the true average abundance of the intracluster 

medium (ICM) is crucial for estimates of its global metal content, which 

in turn is linked to its past evolution and to the star formation history 

^ . of the stellar component of the cluster. We analyze here how abundance 

f~-. | gradients affect commonly adopted estimates of the average abundance, 

assuming various plausible ICM density and temperature profiles. We 
find that, adopting the observed abundance gradients, the true average 
mass weighted abundance is less than (although not largely deviating 
from) the commonly used emission weighted abundance. 
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Q ■ 1. Introduction 

The metal content of the ICM has proved to be a powerful tool to constrain 
the past supernova history of the stellar population of galaxy clusters, being it 
directly linked to the total number of supernovae exploded in the cluster (e.g., 
Renzini et al. 1993). Supernovae also provide a heating mechanism of the ICM, 
via galactic winds that they power, and the relevance of this heating is a much 
debated topic in recent studies of ICM evolution (e.g., Ponman, Cannon, & 
Navarro 1999; Wu, Fabian, & Nulsen 2000; Lloyd-Davies, Ponman, & Cannon 
2000). So, estimates of the mass of metals in the ICM (M^ M ) are precious pieces 
of information, and the more accurate they are, the more constraining they can 
be in recovering the past supernova activity and from here the past histories of 
star formation and heating and evolution of the ICM. A quantity strictly related 
to M^ CM is the average abundance <Z\qm>= M^ cm /Miqm- There are two main 
possibilities to estimate <Z\qm>- One is given by: 

_ /9lCM(x)ZiCM(x)(i 3 X 

I 1 I <^ICM>dopr : 

P1Cm(x) d 6 X. 

where Picm( x ) and Zicm(x) are three-dimensional gas density and abundance 
distributions (e.g., derived from deprojection of observed two-dimensional quan- 
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tities). The other is an emission weighted average abundance: 

r n? CM (x) A[T(x), Z(x)] Zicm(x) d 3 x 



<-^ICM>J 



^ CM (x)A[T(x),Z(x)]d 3 x 



where A is the cooling function. This <Z\qm>l is estimated for example as a 
fit parameter when modeling the observed X-ray spectrum of the whole ICM. 



2. The problem 

There are a few serious problems with estimating <Zicm> using method |lj. In 
fact, for many/most clusters we do not know: 

- the shape of the ICM distribution; 

- the viewing angles under which we are observing the ICM distribution. 
So, we cannot unambiguously deproject the observed quantities to derive 3-D 
ones such as Zicm( x )> Picm( x ) and Ticm( x ). In addition, it is well known that 
deprojection is a very demanding process, sensitive also to: 

- the properties of the instrumental PSF; 

- the measurement errors (e.g., Finoguenov & Ponman 1999). 



Correspondingly, there are clear advantages with method 2 

- the knowledge of the true shape of the ICM distribution is not required 
and the result is independent of the viewing angles; 

- it does not suffer from instrumental effects and amplification of measure- 
ment errors. 

- the average <Z\cm>l is the only accessible information that we may have 
for distant clusters or those for which there are not counts enough for a spatially 
resolved spectroscopy (many more of these are likely to be observed/discovered 
soon, thanks to Chandra and XMM). 



Unfortunately, method 2 has one problem, whose relevance has not been inves- 
tigated so far. Emission weighted average abundances <Zicm>l derived from 
observed X-ray spectra are equal to the true <Z\qm> only in case of a spatially 
independent metal distribution. So, only if Z(x) = Zq constant, then <Z\qm>l 
=<Zicm>= Zq. Recently spatial abundance gradients have been derived for 
the ICM of many clusters from ASCA and BeppoSAX data (e.g., Finoguenov, 
David & Ponman 2000; Irwin & Bregman 2001; De Grandi & Molendi 2001). 
When such gradients are present, in general < Z\qm> L ^ < Z\cm> ■> in a way 
dependent on the spatial distribution of picm, ?icm, ^icm- Here we address the 
following point: how much discrepant are <Z\qm>l and <Z\cm> ? 

Our approach to find the answer is to "calibrate" the <Ziqm>l / <Zicm> 
ratio using many plausible (spherically symmetric) profiles for picm, ^icm and 
Zicm- The ingredients entering the estimates of <Z\cm>l and <Z\cm> are the 
following: 

1) the cooling function A = A(T, Z) over (0.5-10) keV; this has been calculated 
with the thermal X-ray emission code of J. Raymond (that gives the same results 
as the MEKAL model within XSPEC for kT > 3 keV). 
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Figure 1. Abundance profiles adopted in this work 



2) An abundance profile described by ^icm(^) = Zq x [1 — (r/r v ; r ) 6 ], where b is 
a free parameter (Fig. 1) and r v [ r is the cluster virial radius. 

3) Various picm and Ticm profiles (Figs. 2 and 3). In a first set of models we 
assume the ICM to be in (non self-gravitating) hydrostatic equilibrium within 
a chosen cluster (dark matter) potential well, both in the isothermal and poly- 
tropic case. In a second set of models we consider a cooling flow description. 

3. Models of ICM in hydrostatic equilibrium 

3.1. Description 

In this first set of ICM models we explored: 

a. The "NFW" cluster potential well (Navarro, Frenk & White 1996), 
where the density profile of the gravitating matter is of the form p c \(r) = 
4/? c i (r s )rf/r(r + r s ) 2 . The isothermal and polytropic ICM density profiles are 
given by: 



PicmO) =Poxexp 



ump 
kTa 



pol 



> Picm^Hpox 



7 kT o 



i 

7-1 



In the polytropic cases the temperature profile is obviously given by 



r pol 
! ICM 



(r) = T x 



Po 



b. The power law models, in which the cluster density profile is given by 
Pci(f) = Pn x (r/r n )~ a , with a < 2. We have now 



Picm(0 = Po x exp 



fim p 
~kT n ' 



pol 



<t> ), PlCM( r )=/ ? x 



1 " ( — f^ 



7-1 



c. The "standard" /3-models (e.g., Cavaliere & Fusco Femiano 1976), where 
the surface brightness profile of the hot ICM is described by Ex(-R) = ^xo x 



l + {R/Rc 



0.5-3/3 



and (3 = 0.5 - 0.8 (Mohr et al. 1999, Jones & Forman 



1999). The (untruncated) ICM density profiles are given by: 



PicmM = Po x 



R c 



-3/3/2 



pol 



PicmM = Po x 



R c 



i f(fi,i) 
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Figure 2. The density profiles of hydrostatic equilibrium models de- 
scribed in section 3, in the isothermal case. 





L °8 Pioi( r )//W°) 




T IC w(r)A ICM (0) 


u 

0.5 
-1 

1.5 
-2 

?5 


[^riJ- | ■ L J_ I ' ' 

rcusp(NFW)^\ 
: fl-models 

r power law 

F- i i i I i i i I i i 


1 1 ' ' ' 1 ' ' '- 

rTw 


1 

0.8 
0.6 
0.4 
0.2 


- ^X^O^S^ ^ N a=0.5 - 

: cusp (NFWp^ssa^^-x : 

- fl-models "° =:::: ^~^-^ :5 --_: 
: power law **/^^^T~ ,; >~ : 

I* i i i I i i i I i i i I i i i I i r>^ 



0.2 0.4 0.6 0.8 
r / r »i 



0.2 0.4 0.6 0.8 
r / r vi 



Figure 3. The density and temperature profiles of hydrostatic equi- 
librium models described in section 3, in the polytropic case. 



where /(/?, 7) can be found in Ettori (2000). Both density distributions are then 
truncated at a radius Rt = 5R C , with 8 a free parameter. 

Constraints on the parameter values have been derived from observational 
results and from the relations M200 ~ < Tx > (Navarro et al. 1996), < Tx > - 
r 200 (Evrard et al. 1996) and a 2 - < T x > (Xue & Wu 2000). 



3.2. Results 

We estimated <Ziqm>l / <Z\cm> by adopting the following parameter values: 
1) for the abundance profile b = 0.2, i.e., we considered a sort of "cuspy" 
profile, as observed in some clusters (e.g., De Grandi & Molendi 2001; Arnaud, 
this meeting), and b = 1, a linear dependence on radius. The results are almost 
independent of Zq. 
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2) the temperatures are Tq = 3 and 10 keV in the isothermal case; these 
are central temperatures, in the polytropic case. This assumption encompasses 
the values observed for rich galaxy clusters (e.g., Markevitch et al. 1998). 

3)l<7<5/3in the polytropic case. This includes the 7 = 1.2 value that 
gives a "good" fit of some observed temperature profiles for ICM's described by 
a (3- model (Markevitch et al. 1998). 

4) for the power law models a = 0.5 and 1.5. 

5) for the /3-models: (3 = 0.5 and 0.8, 5 = 5 and 10. 

The general trend of the values derived for <Zicm>l / <Zicm> is that: 

• <Z\cm>l I <-^icm> always increases for centrally steeper abundance profiles 
(i.e., lower b values) and steeper gas density profiles (i.e., higher To, a, (3 and 5). 

• <Zicm>l I <Z\qm> is always > 1, and it assumes the values below: 

isothermal case polytropic case 

a. NFW: 1.4-1.7 1.4-1.7 

b. POWER LAW: 1.3-1.8 1.3-1.5 

c. /3-MODELS: 1.4-2.1 1.5-2.2 

• The typical range of values for <Z\qm>l / <Z\cu> is 1.4-2, quite insensitive 
to variations of density and temperature profiles in the chosen (large) range. 



4. Cooling flow models 

We assume the ICM to be described by two gas phases, at two fixed temperatures 
T coo i and T amD , in pressure equilibrium within the cooling radius r coo i. Outside 
r coo i there is only the ambient gas at T am b (Ettori 2001). So: 



PICM(^) =P0,cool x 



/ r 



\ '"cool 



l-5/3 coo i 



+ P0,amb X 



, r 

1 + 1- 
7V 



-1.5,9a 



Correspondingly the X-ray surface brightness profile is given by the superposi- 
tion 



^x(-R) = Sxo,cool x 



( R 



\ '"cool 



0.5+3/3 cool 



+ £ 



X0,amb 



1 + |^' 2 

r r . 



0.5-3/3 amb 



For two well studied cooling flow clusters, for which this decomposition has 
been made (Ettori 2001), we find <Ziqm>l / <Z\cm> as given in Table 1 (for 
b = 0.2 and 1). 



5. Future developements 

We plan to do the following further investigation: 

1) Calibrate <Z\qm>l / <Z\cm> f° r axisymmetric and triaxial models (not 
significantly different results are expected). 
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Table 1. Results for cooling flow models 



Cluster 


O T c — '"cool 

(/i^Mpc) 


J cool 

(keV) 


-'amb 

(keV) 


Pcool 


Pamh 


<^ICM>L / <^ICM> 


A1795 
A2199 


5.8 0.26 
11.5 0.13 


0.5-1.1 

0.4-0.8 


7.4 
4.6 


1.801 
1.635 


0.761 
0.644 


1.6-2.0 

1.7-2.1 



2) Calibrate <^icM>depr / <-^icm> for axisymmetric and triaxial models, by 
considering their projection at different viewing angles, by circularizing and de- 
projecting their X-ray properties, by deriving putative 3-D picm(^) and ^icm(^) 
and comparing M^, Ud = Air J r 2 picM(r)ZicM(r)dr with the true value Mi CM . 

3) Evaluate < Zjqm >l / < -^icm > for ICM's resulting from high resolution 
numerical hydrodynamic simulations that include dark matter, gas and star 
formation (Cen &: Ostriker 2000). 
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